Upper limit on the critical strength of central potentials in relativistic quantum 

mechanics 



Fabian BratQ 

Groupe de Physique Nucleaire Theorique, Academie Universitaire Wallonie-Bruxelles, 
Universite de Mons-Hainaut, B-7000 Mons, Belgique 

In the context of relativistic quantum mechanics, where the Schrodinger equation is replaced by 
the spinless Salpeter equation, we show how to construct a large class of upper limits on the critical 
value, gi , of the coupling constant, g, of the central potential, V(r) = —gv(r). This critical value 
is the value of g for which a first £-wave bound state appears. 
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O I. INTRODUCTION 

(N 

O I A covariant description of bound states of two particles is achieved with the Bethe-Salpeter equation This 
equation reduces to the spinless Salpeter equation j3| when the following approximations are performed: 

• elimination of any dependences on timelike variables (which leads to the Salpeter equation Q). 

• any references to the spin degrees of freedom of particles are neglected as well as negative energy solutions. 
The spinless Salpeter equation takes the form (h = c = 1) 

(N 
O 



o 
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vV + to 2 + V(r) *(r)=M*(r), (1) 



where m is the mass of the particle and M is the mass of the eigenstate (M = m + E, E is the binding energy). We 
restrict our attention to interactions which are introduced in the free equation through the substitution M — > M— V(r), 
where V(r) is the time component of a relativistic four-vector. The interaction could also, in principle, be introduced 
through the substitution p — > p — A(r), where A(r) is the spatial component of a relativistic four- vector. However we 
do not consider this kind of potentials since the derivation of the spinless Salpeter equation from the Bethe-Salpeter 
equation leads to A(r) = 0. Equation JQ) is generally used when kinetic relativistic effects cannot be neglected and 
when the particles under consideration are bosons or when the spin of the particles is neglected or is only taken 
into account via spin-dependent interactions. Despite its apparent complexity, this equation is often preferred to the 

S Klein-Gordon equation. Equation appears, for example, in mesons and baryons spectroscopy in the context of 
potential models (see for example [j, 0, lj, LA M, ) [For a review of several aspects of the "semirelativistic" description 
!> . of bound states with the spinless Salpeter equation see: W. Lucha and F. F. Schoberl, Int. J. Mod. Phys. A 14, 2309 
(1999) and references therein.] 

Due to the pseudo-differential nature of the kinetic energy operator, few exact results are known about this equation. 
■ Most of these results have been obtained for a Coulomb potential (for example, upper and lower bounds on energy 
levels) HE [H[ll Il 3| . Recentl y, u pper and lower limits on energy levels have also been obtained for some other 
particular interactions [lj, [TjJ Hs 03 • 

Conversely to the Schrodinger equation, for which a fairly large number of results giv ing both u p per and lower limits 
on the number of bound states can be found in the literature (see for example [TsL llH. Efl l2ll 122 . l23l Izi . |2"fiL l2?l| ) . 
only two results are known for the spinless Salpeter equation |28[ l2flj. The first result, obtained in Ref. (28j|T is an 



upper bound on the total number of bound states yielding a lower limit on the critical value, gi°\ of the coupling 
constant (strength), g, for which a first S-wave (£ = 0) bound state appears (£ being obviously the angular momentum) 
in the potential V(r) = —gv(r). The second results, obtained in Ref. |29j, is an upper limit on the number of ^-wave 
bound states which yields a lower limit on the critical value, g^ for which a first £-wave bound state appears. 

In this article, we obtain accurate upper limits on the critical strength g^ applicable to attractive (purely negative) 
central potentials which are less singular than — at the origin. This limitation has a deep reason. Indeed, it is 
known that for the spinless Salpeter equation, a potential which behaves like — r _1 at the origin is characterized by a 
maximal value of the coupling constant above which the spectrum is no longer bounded from below. This particularity 
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has been studied in detail for the Coulomb potential (see for example j9|). The — r^ 1 singularity is a critical singularity 
for the spinless Salpeter equation just as the —r~ 2 singularity is a critical singularity for the Schrodinger equation. 
So in this article we discard this class of potentials which should be treated separately. Moreover we suppose that 
the central potential V(r) is piecewise continuous for r G ]0,oo[. The upper limits on we obtain in sec. ITT1 arc 
compared with the exact critical value obtained numerically for some test potentials. These comparisons indicate that 
the new upper limits are very restrictive. Some conclusions are presented in sec. 11111 



II. UPPER LIMIT ON THE CRITICAL STRENGTH 



The idea used to derived the upper limit on is to transform the standard eigenvalue problem obtained with the 
time independent spinless Salpeter equation Q), and where the eigenvalues are the eigenenergies, into an eigenvalue 
problem where the eigenvalues are the critical coupling constants. These critical values of the strength of the potential 
correspond to the occurrence of an eigenstate with a vanishing binding energy. We thus consider the zero binding 
energy spinless Salpeter equation that we need to write as an integral equation. This has be done in Ref. p9l | but 
since we need some modifications in the development, we recall the main line here. 

We have to calculate the Green function of the kinetic energy operator. Similar calculations have also already been 
performed previously pffll l3ll | . In contrast to results found in Ref. [HI , we need here to calculate the Green function 
of the following operator 



T (p 2 ) = v/p 2 + to 2 - m. 



This is done by performing the integral 

G(m,A) 

where A = r — r' and A = |A|. We find that 



(27T) 



dp 



exp(— i p • A) 
\/ p 2 + m 2 — m 



(2) 



(3) 



G(m, A) = 



4ttA 



1 + -F(toA) 
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F(toA), 
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JL Kl {z) + - = K l {y) + -- J dzK (z), 



(4a) 



(4b) 



and where K u (y) is a modified Bessel function (see for example [321 p. 374]). The zero binding energy spinless 
Salpeter equation takes thus the form of the following integral equation 



*(r) 



dr'G(m,A)F(r')#(r'), 



(5) 



with G(m, A) given by Q. We now restrict our attention to central potentials V(r) = V(r), with r = |r|. 

Integration over angular variables reduces the integral equation JSJ to the following one-dimensional integral equa- 
tion 



with 



u £ (r) 



dr' G e (m,r, r')V{r')u e (r'), 



d0> sm9' H( -™ A ^ Pi(cosfl'), 



(6a) 



(6b) 



where ue(r) is the radial wave function, \&(r) = (ui(r)/r)Yg m (r ) and where H(x) is defined by (|4a|l . 

An important technical difficulty, to obtain a symmetrical kernel, appears if the potential possesses some change 
of sign (see relation Q below). This is overcome when one searches for necessary conditions, or upper bound on the 
number of bound states, by replacing the potential by its negative part V(r) — * V~(r) = max(0, —V(r)). Indeed, the 
potential V~(r) is more attractive than V{r) and thus a necessary condition for existence of bound states in V~(r) is 
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certainly a valid necessary condition for V(r). This procedure can no longer be used to obtain sufficient conditions. 
For this reason we consider potentials that are nowhere positive, V(r) — —gv(r), with v(r) > 0. 

The integral equation JHJ can be written with a symmetrical kernel provided we introduce a new wave function 



Mr) = \V{r)\ 1 ' 2 u t {r). 
This change of function leads to the following integral equation 

/>oc 

<M r ) = 9 dr' Kt(m,r,r')<t>t(r'), 
Jo 

with 

K e (m, r, r 1 ) = w(r) 1 / 2 G t (m, r, r') w(r') 1/2 . 



(7) 



(8b) 



The relation JHJ) is thus an eigenvalue problem and, for each value of I, the smallest characteristic number is just the 
critical value g c • The other characteristic numbers correspond to the critical values of the strength of the potential 



for which a second, a third, ^?-wave bound state appears. The kernel (|8b|) acting on the Hilbert space L 2 (K) is 
an Hilbert-Schmidt operator for potentials which decrease faster than r _1 at infinity. Thus this kernel satisfies the 
inequality 



dxdy K e (x,y)Ke(x,y) < 



(9) 

pp. 102-106] (in 



Consequently the eigenvalue problem JHJ always possesses at least one characteristic number 
general, this problem has an infinity of characteristic numbers) . 

Now we use the theorem (see for example (3^, PP- 118-119]) which states that, for a symmetric Hilbert-Schmidt 
kernel, we have the variational principle 



max 



for tp(r) satisfying 



dx dy K e (x,y) (p(x)ip(y) 



dr tp(r) 2 = 1. 



1 



(10) 



(ii) 



The maximal value is reached for (p(x) = ip\{x), where ipi(x) is the eigenfunction associated to the smallest eigenvalue 
g\. Consequently, for an arbitrary normalized function, f(x), we obtain the following upper limit on 51 



dxdyK e (x,y)f(x)f(y) 



(12) 



For the clarity of the discussion we now consider in two separate sections the ultrarelativistic regime (m = 0) and 
the relativistic regime (m > 0). 



A. Ultrarelativistic regime m = 

In this section, we derive an (among others) upper limit on the critical value, gip , of the coupling constant, g, of 
the potential, V(r) — —gv{r), for which a first £-wave bound state appears in the ultrarelativistic regime (m = 0). In 
this limit, the kernel takes a simple form since mK\{my) = X/y when m goes to zero. This implies that 

2 

lim mHtmA) = —. (13) 

m^O 7rA 

The function Gt(0, r, r') takes then the form 

G i (0,r,r') = ^- f d9'^P e (cose'). (14) 



7T 
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A simple change of variable leads to (3^, P- 335] 

Ge(0,r, r') = i- £ 



(r 2 + r' 2 )/{2rr') - y 



1 

-( 

n 



2rr' 



(15) 



where the function Qt{x) is a Legendre function of the second kind. The function Gi(0,r,r') can thus be evaluated 
explicitly for each value of the angular momentum t. We have for example 



Go(0,; 



= Iln 

7T 



r + r' 



and 



Gi(0,r,r') 



,.'2 



2rr' 



In 



(16) 



(17) 



Since the function Gg(0,r,r') is given by the relation l|15|) . it follows that the kernel -Kf(0, r, r'), see 1)8 bfl . is known 
for each value off. Now, we just need to choose a suitable normalized function f(r) to apply the variational principle 
reported above. 

For simplicity we restrict the rest of the following discussion to I — but extensions to non vanishing values of the 
angular momentum is obvious, one just need to computed the corresponding expression of the function Gf(0, r, r'). 

The function f(r) should be as close as possible to the zero binding energy wave function but also should be general 
and simple enough to obtain a neat formula. We simply choose 



/(r)=A[r p - 1 v(rf] 1/2 , p > 0, 



(18) 



where A is the normalization factor. The relations JSEJ, i|12|) . I|l(j|) and 118|) lead to the following upper limit on g, 



(0) 



7 (°) < 



a7r L°° dx Fi(2p — l;x) 



2 Fx (p; z) dy F, (p; y) In (f±f ) 



= r7 m=0 
— i/up,l ' 



(19) 



where Fi(g;x) = x^ 9_1 ^ 2 u(x)^ +1 ^ 2 and where we have introduced the parameter a which takes the value 1 respec- 
tively 2 for one respectively two (identical) particle problems. The most stringent upper limit is obviously obtained 
by minimizing the right hand side of l|19|) with respect to all positive values of p. 

A simpler, but less stringent, version of this upper limit can be obtained with the help of the following minorization 



In 



> 



2y 



(20) 



,(o) 



< 



cot J °° dx F\ (2p — 1; x) 

4 Jo°° dx x ~ 1 Fi (p> x ) Jo d v y Fi (p> y) 



5up,2 



(21) 



The accuracy of these upper limits can be tested with some typical potentials. The comparison between the exact 
results (obtained by solving numerically the spinless Salpeter equation) and the upper limits (|19(l and (|21|l is reported 

in Table [I] We have also added two lower limits on gi°^ obtained with the upper limits on the number of bound states 
derived in Refs. [28ll2^| . Note that for these tests, we choose a two identical particles problem, a = 2. 

The results reported in Table [I] indicate clearly that the accuracy of the upper limit l|19l) is quite remarkable. The 
upper limit l|21|) is obviously less stringent but could prove to be useful to obtain explicit formulas. The typical value 
of p which optimize these upper limits varies between 2 and 3. We do not consider other choices for f(r) (see JTS}) 
here since the relation l|19|) is already very accurate. 

As an additional indication that the upper limits obtained with the method proposed in this work are quite accurate, 
we report in TablclTTIa comparison between the exact value of the critical strength gi 1 ^ {i = 1) and the corresponding 
upper limit obtained with the relations |Q , ifT7|l and (|T%)l and noted 3^ =0,£=1 in this Table. 
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B. Relativistic regime m > 

(£) 

To obtain an upper limit on <?c for a non vanishing mass m, we need to calculate the expression of the function 



Gi(m, r,r'). To this end we note that 



K 1 (y)<F(y)<K 1 (y) + ^, (22) 



^Frorn the relation l|12l) . it is obvious that a minorization of the kernel Ki{m,r,r') is enough to obtain the upper 
limit. However, the minorization (|22|l of the function F(y) is too crude to obtain good results. Instead we use 

F(y)>K 1 {y) + ~-~ ex p{-y). (23) 

This minorization l|23[) is proved in Appendix 1X1 From the definition of Ge(m,r,r') l)6b|) and the inequality l|23|) we 
obtain 

Ge(m 7 r, r') > —Qi{m, r, r') + Se(m, r, r') - -%(m, r, r'), (24) 

7T Z 

where 

/r+r' / 2 i 12 2 \ 

f dyK 1 (my)P l \l ^~ V j , (25) 

t-r+r' / r 2 +r i2_ y 2 

Si{m,r,r) = ml dy P t 



with r< = min[r, r'] and r> = max[r, r'] and 



r - r i\ \ 2rr' 



2m 
2£+ 1 



r*< +1 r><, (26) 



i exp 



(—my/ r 2 + r' 2 — 2rr'yj 



%>(m,r,r') = mrr' / dy , „ L Pe(y), 

yr + r — zrr'y 



i X"i 



a/2 ( —my/ r 2 + r' 2 — Irr'y ) 
= \lzm 2 rr'l dy ^ - ^ _ P,( y ), 

= 2m\fr? K i+ 1 (mr>) 7g + 1 (mr<), (27) 

where /^(x) is a modified Bessel function (see for example fH, p. 374]). The kernel Sg(m,r,r') is actually the Green 
function of the nonrelativistic kinetic energy operator and takes a simple form while the kernel G(,(m,r,r') can be 
calculated analytically for each value of £ [3(J, |31j • We find for example 

So(m,r,r') = K (m\r - r'\) - K (m(r + r')), (28) 
C?i(to, r, r') = i^o(w|r — r'|) + Ko{m{r + r')) 

+ [(r + r') Ki(m(r + r')) - \r - r'\ KUmlr - r'\)] . (29) 

mrr' 

Now, we just need to choose a suitable normalized function f(r) to apply the variational principle reported above 
(see We take the following expression for f(r) 

f(r) = A [r 2p - 1 u(r) p ] 1/2 , p > 0. (30) 

For simplicity we again restrict the rest of the following discussion to £ — but extensions to non vanishing values of 
the angular momentum is obvious, one just need to computed the corresponding expression of the function Gi(m, r, r'). 
The relations (Sb} , (JT^J , (fT5|) and lead to the following upper limit on gi°' 

g{0 ) < _ af™dxy/ZF 2 {2p-l;x) ^ ^ m>Q ^ 

~ 2 fo° dx F 2(p;x) Jo dy F 2 (p,y)T(x,y) up 
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with 

T(x,y) = ~[K (f3(x - y)) - K (f3(x + y))] + 2f3y + ~[exp(-/?(x + y)) - exp(-/?(z - y))], (31b) 

where F2(q,x) = x q ^ 1 ' 2 v(x)^ q+1 ^' 2 and where [3 — mR, R being the scale of length which appears in the potential 
(v(r) — v(Rx)). Again, we have introduced in l|31a|) the parameter a which takes the value 1 respectively 2 for one 
respectively two (identical) particle problems. 

The accuracy of this upper limit can be tested with some typical potentials. The comparison between the exact 
results and the upper limit 1)3 l|l is reported in Table ITTll Note that for these tests, we also choose a two identical 
particles problem, a — 2. 

The results reported in Table II I II indicate clearly that the accuracy of the upper limit l|31|l is quite good. But 
for small value of (3 the upper limit is however less restrictive. Thus for small value of (3 it is preferable to use an 
intermediate form for f(r). We then propose in general to use 

/(r) = A [r ap_1 v(r) p ] ^ 2 , p > 0, (32) 

with 1 < a < 2. This last expression for f(r) improves significantly the restriction on the possible values of 
Indeed, for (3 = 0.1, a — 1.18 and for the exponential potential, the upper limit is then equal to 4.812 instead of 5.390. 
But this additional flexibility is only significant for small value of f3, indeed for (3 = 0.5 the best upper limit is found 
to be equal to 2.521 (for the exponential potential and a = 1.69) instead of 2.547. 

However, even with the choice ()32JI for /(r), the upper limit l|31(l still yields less restrictive results for small (3 than 
those obtained for larger values of (3 or those obtained with the upper limit (|T§|l . This is easy to understand, since 
this is in the sector of small (3 that the error introduced by the inequality (|23|l is the most important. Indeed, in the 
limit of (3 going to zero, the upper limit 1)3 l|l coincides with the upper limit (|19|l and for (3 going to infinity, only the 
non relativistic kernel Si{m, r, r') contributes. 



III. CONCLUSIONS 

In this article we have shown how to construct upper limits on the critical value, g c , of the coupling constant, 
g, of a central potential, V(r) — —gv(r). The method used to derive the upper limits is quite general and other 

(possibly more complicated) families of upper limits yielding (possibly) stronger restrictions on could also be 
obtained. Indeed, the method is based on a variational principle for which a trial zero energy wave function is needed. 
There is no limitation on the accuracy of such a trial function, which imply that there is, in principle, no limitation 

(€)... 

on the accuracy of the upper limit on g c derived with this procedure. However, this remark is only true for the 
ultrarelativistic regime, m — 0, where the kernel of the integral equation has been calculated exactly. For m > 0, a 
minorization of the kernel has been used yielding some errors in the restrictions on the possible values of the critical 

(£) 

value g^ which cannot be compensated by a better choice of the trial zero energy wave function. In this article we 
have proposed in sec. [H] a compromise between accuracy and simplicity of the final formula. The accuracy of the 
upper limits on gi^ was then tested with some typical potentials. 
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APPENDIX A: MAJORIZATION OF THE PRIMITIVE OF K (x) 
We choose the following integral representation for the modified Bessel function Kq(x) [22, p. 376] 

/>oo 

K (x) = / dt cxp(-xcosht). (Al) 
Jo 

We have 

f°° , „ , s f°° , exp(-ycoshi) . , f°° , 1 tt . , 
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TABLE I: Comparison, for some typical potentials, between the exact critical values, g^' , the upper limits <?u^>,i U 1191 . <?2p 
(121 1 and the lower limits obtained in Refs. |28U29|. 



v(x) 


Ref. [29] Rcf. [28 


1 


m=0 
Sup.l 


m=0 
Sup, 2 


cxp(-x) 


4.443 


4.370 


5.574 


5.584 


7.411 


[cosh(x)]~ 2 


4.126 


3.886 


5.008 


5.018 


6.769 


exp(— x 2 ) 


4.513 


4.169 


5.426 


5.442 


7.399 


x exp(— x) 


3.696 


3.349 


4.360 


4.364 


5.964 



TABLE II: Comparison, for some typical potentials, between the exact critical values, g c and the upper limit g, 
obtained with the relations 112H . 1171 and 1181 . 



v(x) 


o (1) 


m=o,e=i 

i/up 


exp(— x) 


10.975 


10.992 


[cosh(s)] -2 


8.1174 


8.1268 


cxp(— x 2 ) 


10.200 


10.231 


x exp(— x) 


9.5442 


9.5636 



TABLE III: Comparison, for some typical potentials, between the exact critical values, gi and the upper limit g^° 13H 



v(x) 


exp( 


-x) 


[cosh (a;)] 2 


exp( 


-x 2 ) 


p 




yu m P >0 


(o) 
9 c 


„m>0 

yup 


o (0) 
yc 


m>0 

yup 


0.1 


4.694 


5.390 


4.461 


4.994 


4.927 


5.363 


0.5 


2.387 


2.547 


2.766 


3.006 


3.309 


3.589 


1 


1.361 


1.407 


1.742 


1.843 


2.198 


2.352 


2 


0.7133 


0.7206 


0.9598 


0.9862 


1.257 


1.307 


3 


0.4804 


0.4817 0.6549 


0.6642 


0.8669 


0.8880 


4 


0.3607 


0.3615 


0.4956 


0.4994 0.6589 


0.6694 


5 


0.2890 


0.2893 


0.3981 


0.3999 


0.5305 


0.5364 



